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We study Gamow- Teller strength distributions in the neutron-deficient even isotopes 184 ~ 194 p D in 
a search for signatures of deformation. The microscopic formalism used is based on a deformed quasi- 
particle random phase approximation (QRPA) approach, which involves a self-consistent quasipar- 
ticle deformed Skyrme Hartree-Fock (HF) basis and residual spin-isospin forces in both the particle- 
hole and particle-particle channels. By analyzing the sensitivity of the Gamow- Teller strength 
distributions to the various ingredients in the formalism, we conclude that the /3-decay of these 
isotopes could be a useful tool to look for fingerprints of nuclear deformation. 
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I. INTRODUCTION 



Neutron-deficient Pb isotopes turn out to be a unique laboratory to explore the phenomenon of shape coexistence 
in nuclei and have recently been a subject of much experimental and theoretical interest. The existence of at least 
one low-lying excited + state in all even-even Pb isotopes between A — 184 and A = 194 has been experimentally 
observed at excitation energies below 1 MeV (see Ref. [1] and references therein). The most extreme cases are 186 Pb 
and 188 Pb with two excited + states below 700 keV [2,3]. These are good examples where the ground state and 
the first two excited states, which are below 1 MeV, correspond to macroscopical shapes. Experimental data on the 
t— I | excited states of these Pb isotopes are not just limited to a few + states. Rotational bands built on top of these + 
states have also been observed [1,4], whose properties have been used to interpret them as corresponding to oblate 
| and prolate deformations. 

The existence of the + low-lying states arises from the combined effect of a proton shell gap at Z = 82 and a 
large number of hole neutrons in the N — 126 core. The influence of the magic proton number Z = 82 is particularly 
strong and the ground state of Pb isotopes is known to be spherical down to 182 Pb. The weakening of the magicity 
of the Z = 82 shell manifests itself through the appearance of low-lying + states. 

Different types of models have been invoked to explain the coexistence of several + states at low energies. In a 
shell model picture, the excited + states are interpreted as two-quasiparticle and four-quasiparticle configurations [5] . 
Protons and neutrons outside the inert core interact through pairing and quadrupole interactions to generate deformed 
structures. Within a mean field approach, this is interpreted as a macroscopic change from a spherical to a deformed 
shape. The energy of the different shape configurations can be calculated using a nuclear potential, where the energies 
of the single particle orbitals depend on the deformation. Potential energy surface calculations have become more 
and more sophisticated with time, resulting in accurate descriptions of the nuclear shapes and the configurations 
involved. Calculations based on phenomenological mean fields and Strutinsky method [6,7], predict the existence of 
several competing minima in the deformation energy surface of neutron-deficient Pb isotopes. Self-consistent mean 
field calculations [8,9] and calculations including correlations beyond mean field [10-12] confirm these results. 

Sophisticated configuration mixing calculations [10-12] have been shown to support the mean field result that these 
neutron deficient systems display more than one minimum corresponding to spherical and deformed shapes. The 
lowest excited states in configuration mixing calculations have average deformations close to the deformation of the 
mean field minima [10-12]. All the approaches analyzed in Ref. [11] from mean field up to very sophisticated angular 
momentum projected generator coordinate method with Gogny forces provide the same underlying basic picture 
of strong coexisting spherical, oblate and prolate shapes. Although the mean field approach does not reproduce 
quantitatively the experimental relative energies of the + states, it provides a good qualitative description of the 
three minima and is a good approximation to associate the + states observed at low energies with coexisting energy 
minima in the energy surface. The physical + observed states have been also considered [13] as superpositions of 
spherical, oblate or prolate configurations, with relative weights determined phenomenologically, in the spirit of the 
shape mixing picture [13]. Such kind of calculations have been done in Ref. [14] for 190 - 200 pb isotopes. 
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Mean field predictions of deformed ground states in the neutron-deficient Pb region are very sensitive to fine details 
in the calculations, specially to pairing effects. This has been shown for example within deformed relativistic mean 
field (RMF) calculations. Standard RMF calculations [9,15] do not reproduce the spherical ground states in this mass 
region, being at variance with experimental data. In RMF, the ground states of neutron-deficient Pb isotopes are 
found to be deformed for standard forces with constant pairing gaps [16]. However, it has been shown [17] that using 
constant strengths for the pairing interaction (which makes the gap parameters strongly dependent on deformation), 
produces spherical ground states for all Pb isotopes. More recently [9], an improved pairing treatment has been used 
by means of relativistic Hartree-Bogoliubov (RHB) calculations, using the NL3 effective interaction and a finite range 
Gogny interaction to describe the pairing properties. The results of these calculations show that 184 Pb and 186 Pb have 
spherical ground states with low-lying oblate and prolate minima in qualitative agreement with experiment. However, 
with increasing neutron number, the oblate minimum is lowered in energy and 188 - 194 Pb have oblate ground states, 
contradicting the experimental data [1]. In order to recover the spherical ground states, a new parametrization of the 
effective interaction was proposed, which takes into account the sizes of spherical shell gaps. 

The special sensitivity of the ground states in this mass region to deformation and pairing effects has been also 
studied from non-relativistic mean field calculations. In Ref. [18] it was found that the potential energy curves in 
neutron-deficient Pb isotopes are very sensitive to the relative intensity of the neutron versus proton pairing strengths. 
In this reference it was shown that slightly different choices of the pairing strengths produce different ground state 
shapes. 

On the other hand, it has been shown [19,20] that the decay properties of unstable nuclei may depend on the nuclear 
shape of the decaying parent nucleus. In particular, the Gamow- Teller (GT) strength distributions corresponding to 
/3 + -decay of proton-rich nuclei in the mass region A s=a 70 has been studied systematically [20] as a function of the 
deformation, using a deformed HF+BCS+QRPA approach. In Ref. [20] the GT strength distributions of several nuclei 
were identified to be particularly sensitive to the deformation of the /3 + -emitter, and a comparison of the theoretical 
results with recent accurate measurements [21] has been used to determine the nuclear shape in neutron deficient Kr 
and Sr isotopes. 

In this paper we study the GT strength distributions of neutron-deficient Pb isotopes from A = 184 up to A = 194 
and their dependence on deformation. The aim here is to identify possible signatures that allow one to determine the 
existence of spherical, prolate or oblate nuclear shapes, in the spirit of the shape coexistence picture, by looking at the 
/3 + -decay pattern. Although possible effects of the quadrupole shape mixing cannot be discarded and a QRPA-type 
approach might not be a good approximation to treat systems undergoing large amplitude collective motions, the 
calculations performed in this paper are perfectly reasonable as a first estimate of the /3 + -decay patterns. This study 
is timely because the possibility to carry out these measurements is being considered at present [22]. 

The paper is organized as follows. In Sec. II we present a summary of the theoretical framework used in our 
calculation. Sec. Ill contains our results on the Gamow- Teller strength distributions of the Pb isotopes, as well as a 
discussion of their dependence on deformation. The conclusions are written in Sec. IV. 

II. THEORETICAL FRAMEWORK 

In this Section we summarize briefly the theory involved in the microscopic calculations. More details can be found 
in Ref. [20] . The theoretical formalism used here is based on a deformed Hartrce-Fock mean field calculation, using an 
effective two-body Skyrme interaction and including pairing correlations in BCS approximation. The single particle 
energies, wave functions, and occupation probabilities are generated from this mean field. We consider in this paper 
the force Sk3 [23], as an example of one of the most simple, traditional and successful Skyrme interaction and the 
force SG2 [24], which has been successfully tested against spin and isospin excitations in spherical [24] and deformed 
nuclei [20,25]. Some results will also be given for the more recent force SLy4 [26]. 

The formalism is safely restricted to axially deformed nuclear shapes because it has been shown [6,12] that no 
features of particular interest, such as minima or flat regions, appear for triaxial shapes in neutron-deficient Pb 
isotopes. Time reversal and axial symmetry are then assumed here. The single-particle wave functions are expanded 
in terms of the eigenstates of an axially symmetric harmonic oscillator in cylindrical coordinates, using eleven major 
shells. We calculate the energy surfaces as a function of deformation for all the even-even neutron-deficient Pb isotopes 
under study here. For that purpose, we perform constrained HF calculations with a quadrupole constraint [27] and 
minimize the HF energy under the constraint of keeping fixed the nuclear deformation. In order to study the above 
mentioned sensitivity of the energy surfaces to pairing correlations, we compare the results obtained using fixed gap 
parameters for protons and neutrons with the results obtained using constant pairing strengths G n . u . The fixed 
gap parameters are determined phenomcnologically from the odd-even mass differences [28] when this experimental 
information is available, otherwise, we use the expression A = 12A- 1 / 2 MeV. The G„, v strengths are determined from 
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To describe Gamow- Teller excitations we use the QRPA adding to the mean field a separable spin-isospin residual 
interaction, which is expected to be the most relevant interaction for that purpose. This interaction contains two 
parts, particle-hole (ph) and particle-particle (pp) . The ph part is responsible for the position and structure of the 
GT resonance [20,29-31] and it is a common practice to fit the strength of the force to reproduce the energy of the 
resonance [31]. The advantage of using separable forces is that the QRPA energy eigenvalue problem is reduced to find 
the roots of an algebraic equation. The particle-particle interaction is a neutron-proton pairing force in the J™ = l + 
coupling channel. We introduce this interaction in the usual way [30,32,33], that is, in terms of a separable force 
with a coupling constant Kq T , which is fitted to the phenomenology. Since the peak of the GT resonance is almost 
insensitive to the pp force, Kq T is usually adjusted to reproduce the half-lives [30,31]. 

The pnQRPA phonon operator for GT excitations in even-even nuclei is written as 

r ^ = E [ X ^ a t a t Y^a,a v ] , (2.2) 

where ir and v stand for proton and neutron, respectively, a + (a) are quasiparticle creation (annihilation) operators, 
uj k are the RPA excitation energies, and ,Y^f the forward and backward amplitudes, respectively. It satisfies 

rv|o> = 0; t+ k \o) = \oj k ). (2.3) 

Solving the pnQRPA equations [29,32], the GT transition amplitudes in the intrinsic frame of even-even nuclei 
connecting the QRPA ground state of the parent nucleus |0) to one phonon states in the daughter nucleus \uok), are 
given by 

(w K \(3±\Q) = T MZ" , (2.4) 

where /3j£ = ctk^, K = 0, ±1, and 

M - K = E («*» x ™ + *™ Y ™ ) : M + K = E («™*™ + 9nuY^ ) , (2.5) 

7TI/ 7TI/ 

with 



q vv = u u v^- q vv = v v u^ ; Y% = {u \a K \ w) , (2.6) 

where v's are occupation amplitudes (u 2 = 1 — v 2 ). 

Once the intrinsic amplitudes in Eq. (2.4) are calculated, the Gamow- Teller strength B{GT) in the laboratory 
frame for a transition I i K i (0 + 0) — ► IfKf(l + K) can be obtained as 

B ± (GT) = \( T f M f 1^1 J * M *>| 2 = {^/.o (^/ l^ko) 2 + 25 Kf ,i (4> Kf \Pt\ 0o) 2 } , (2.7) 

Mi,Mf,n 

in units of g\/4ir. To obtain this expression we have used the initial and final states in the laboratory frame expressed 
in terms of the intrinsic states \<f>K >, using the Bohr-Mottelson factorization [34]. 

The /3-decay half-life is obtained by summing up all the allowed transition probabilities, weighed with phase space 
factors, up to states in the daughter nucleus with excitation energies lying within the Q-window, 

T r/2 = ^E/( Z ^) B ( GT )' ( 2 - 8 ) 

where / (Z, u) is the Fermi integral [35] and D = 6200 s. We include standard effective factors 

A 2 = [( 9 a/9v) cS } 2 - [0.7 (9A/gv) flc f ■ (2.9) 

In (3 + /EC decay, / (Z, u>) consists of two parts, positron emission and electron capture. In this work we have computed 
them numerically for each value of the energy, as explained in Ref . [35] . 
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III. RESULTS AND DISCUSSION 



In this section we present the results obtained from the above mentioned theoretical formalism for the /3-decay 
patterns of the neutron-deficient 184 > 186 ' 188 > 190 < 192 > 194 pb isotopes. First, we discuss the energy surfaces and shape 
coexistence expected in these isotopes. Then, we show on some examples the sensitivity of the Gamow- Teller strength 
distributions to the various ingredients in the calculation. Finally, we present the results obtained for the Gamow- 
Teller strength distributions with special attention to their dependence on the nuclear shape and discuss their relevance 
as signatures of deformation to be explored experimentally. 



We show in Fig. 1 the HF energies calculated with the Skyrme force Sk3 and pairing correlations treated in the 
fixed gap approach, for all the isotopes considered, as a function of the quadrupole deformation parameter (3. The 
latter is defined as 



in terms of the proton quadrupole moment Q p and the charge r.m.s. radius < r 2 >. For a better comparison, the 
curves are scaled to the energy of their ground states and are shifted by 1,2,3,4,5 MeV for A=186,188,190,192,194, 
respectively. As we can see in the figure, prolate solutions corresponding to values (/3 > 0) are always present from 
A=184 to A=190, and arc indeed the ground state for this choice of force and gap parameters. The prolate minimum 
shows only a shoulder in A=192 and disappears completely in A=194, where only the spherical solution remains. 
Spherical local minima can be seen in all the isotopes except for A=192, the one in A=190 being very shallow. Oblate 
minima (f3 < 0) are also present in all the isotopes except for the one closer to stability A=194. We can also see in 
Fig. 1 that the oblate deformation corresponding to the local oblate minima in all the isotopes is rather stable at 
(3 = —0.2, while the prolate deformation at the minima becomes larger as we remove more and more neutrons, to 
reach (3 — 0.32 in A=184. Mass deformations follow closely charge deformations. 

Our results are in qualitative agreement with experiment in the sense that we find different coexisting shapes in this 
mass region, although they are at variance with experiment [1] in the predicted ground state shape. Experimentally, 
the ground states of all the isotopes considered here are found to be spherical, while we find a spherical shape as the 
ground state in A=194, an oblate shape in A=192 and a prolate shape in the remaining cases considered (A=184-190). 
This discrepancy is not surprising because as we have mentioned in the Introduction, even the more sophisticated 
formalisms fail in some instances [9] to account for the spherical ground state, showing an extreme sensitivity to fine 
details of the two-body interactions. To illustrate this sensitivity, we compare in Fig. 2 on the example of A=184, 
the energy deformation curves obtained when we change the force from Sk3 to SG2 and SLy4, and when we change 
the pairing correlations from fixed gap parameters to constant pairing interaction with G n ( v ) strengths given by Eq. 
(2.1). We can see in Fig. 2 that although we get a prolate ground state in all the cases considered, the calculation 
with a constant pairing strength (dotted line) takes the spherical minimum very close in energy to the ground state. 
Indeed by slightly changing the G- values (the relative values of G^ versus G„), we can get the spherical solution as the 
ground state. An analysis along these lines was made in Ref. [18] with analogous results. It is also worth mentioning 
that the force SLy4 produces three minima which are very close in energy. 

The results shown in Fig. 2 represent what happens when using other forces and other gap parameters or pairing 
strength values. That is, the deformation minima remain at about the same deformation values but their relative 
energies may change considerably. This example is also representative of the situation with respect to other isotopes. 



In this subsection we study the energy distribution of the Gamow- Teller strengths calculated at the equilibrium 
shapes that minimize the energy of the nucleus. First, we settle the values of the coupling strength constants of our 
separable residual ph and pp forces. In this work we use the same value for all the isotopes considered, which is 
fixed by comparison to the experimental B~{GT) strength resonance observed from (p,n) charge exchange reactions 
in 208 Pb at an energy of 19.2 MeV relative to the parent nucleus 208 Pb [36]. In this way, we reproduce the energy of 
the resonance with Xgt = 0-09 (0.07) MeV, when using SG2 (Sk3) forces and these are the values we use here. Once 
the value of the ph spin-isospin residual force has been established, the value of the pp residual force is usually chosen 



A. Deformation energy curves 




(3.1) 



B. Gamow- Teller strength distributions and deformation 
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to reproduce the measured half-life. In the present case, we have calculated the half-lives as a function of the k P q T 
parameter and have found that they are very little dependent on this parameter. This is illustrated in Fig. 3, where 
we can see the results obtained for T1/2 from Eqs. (2.8) and (2.9). The results correspond to the case 190 Pb with an 
experimental half-life T X j 2 = 71 s. The effect of the pp force on the GT strength distributions can be seen in the inset, 
where we have plotted the distributions for the oblate and prolate shapes with k P q T — and 0.02 MeV. From now on, 
we use kPq T = 0.02 MeV in the calculations of the GT strengths. This value is compatible with the parametrization 
of Ref. [31], where it was found k p J t = 0.58/A 7 MeV. 

In Fig. 4 we study the influence of both the Skyrme force and the pairing treatment used on the GT strength 
distributions. We do that on the example of the isotope 184 Pb and perform the calculations for the different shapes 
that minimize the energy of this isotope, which are shown in Fig. 2. The purpose of this figure is to show that the 
different profiles obtained for different deformations remain when the force and/or pairing are changed. This implies 
that the GT profile is characteristic of the shape and depends little on the details of the two-body forces used. Fig. 4, 
as well as all the coming figures containing GT strength distributions, shows in addition to the discrete strengths, a 
continuous distribution resulting from a folding procedure using T = 1 MeV width gaussians on the discrete spectrum. 
The B(GT) strength is plotted versus the excitation energy of the corresponding daughter nucleus. On the left panels 
in Fig. 4 we show the results from the oblate shape, in the middle the results from the spherical shape, and on the 
right panels we show those from the prolate shape. For an easier comparison between the various results, we have 
used in all the panels the same scales of strengths and energies. In the upper row we compare the results from two 
different Skyrme forces Sk3 and SG2, using pairing correlations with fixed gap parameters A T — A v — 0.9 MeV, as 
obtained from the general expression given in Sec. 2. Results from Sk3 and SG2 arc plotted upward and downward, 
respectively. In the lower row we compare the results from Sk3 using fixed gap parameters (plotted upward) and 
constant pairing strengths given by Eq.(2.1) plotted downward. The main point to stress in this figure is the fact 
that the structure of the GT strength profiles is different for each deformation but very similar when we compare the 
calculation of reference (Sk3 with fixed pairing gaps) plotted upward with those plotted downward. This is specially 
true in the low energy region where the oblate shape gives rise to a three bump structure below 5 MeV, the spherical 
shape a single bump structure, and the prolate shape a two peaked structure with a first strong single peak at very 
low energy and a more spread second bump at higher energies. 

The main difference between the calculations from Sk3 and SG2 is a small shift to lower excitation energies in 
the case of SG2. This is better seen in the spherical case where a single peak is present. Actually, in the strict 
spherical limit the fragmentation of the strength observed in the figure would collapse to single excitation energies 
corresponding to the transitions between degenerate spherical orbitals. The fragmentation observed arises from the 
small deformation present in our deformed formalism (/3 = 0.003), which is enough to break the degeneracy of the 
third components of angular momentum. The origin of the shift in the excitation energy observed with different 
forces, can be traced back to the predicted energies for the /19/2 and /111/2 spherical shells. The main contribution to 
the GT strength at low energies comes from the transition connecting the almost fully occupied irhii/ 2 shell with the 
partially unoccupied vhg/2 shell, and the relative position of the two shells (closer for SG2 than for Sk3) determines 
the GT excitation energy. This effect translates into the deformed cases where the non-degenerate orbitals spread 
around the energy of the spherical shell. The difference between the calculations with different pairing treatments 
shown in the lower panels is due to the different occupation probabilities of the orbitals that weigh differently the 
matrix elements connecting the proton and neutron states, as well as to the different quasiparticle energies that finally 
shift the excitation energy of the GT transitions. 

In Fig. 5 we show the results for the B(GT) strength distributions in (g^/An) units, calculated within QRPA based 
on the force Sk3 with fixed pairing gaps and residual separable interactions ph and pp with coupling strengths given 
by Xgt = 0-07 MeV and k P q T = 0.02 MeV, respectively. In the left panels we show the results obtained with spherical 
shapes from 184 Pb up to 192 Pb, the middle panels are for oblate shapes and the right panels are for prolate shapes 
in the same chain of isotopes. In each column we have used the same scale for the GT strength to appreciate better 
the changes as we increase the number of neutrons. The scale on energies is always 7 MeV, enough to include all the 
Qec energies, whose maximum value is Qec — 5-84 MeV in 184 Pb. We can see in this figure the progress of the 
profiles as we change the mass number (number of neutrons). We can see a continuous transition from the lightest 
(A=184) to the heaviest (A=192) isotopes, that can be guided by the dashed lines crossing the figures. Thus, we can 
see that the strength in the spherical case is accumulated in a resonance peaked at around 3 MeV and quite stable 
against the number of neutrons. This strength has its origin in the nhu/ 2 — ► vhg/2 transition. In the case of oblate 
shapes, we observe one peak at low energies around 2 MeV, which stays in all the isotopes although is weaker as we 
approach the stability. There is also a second peak clearly seen in 184 Pb that moves to lower energy with mass and 
merges with the first peak gradually. Finally, we observe a third peak that again moves to lower energy with mass. 
While this third peak is practically the same in all the isotopes, the two first peaks become weaker with increasing 
neutron number. An analysis of the structure of these transitions reveals that while the first two peaks arise from 
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several transitions between negative parity states (mainly from 7rft,n/ 2 — * ^9/2 transitions), the third peak at higher 
energy corresponds to transitions between positive parity states with a dominant transition 7r(13/2) + — > v(ll/2) + , 
connecting the proton i 13 / 2 shell with the neutron in/2 shell. The prolate cases show also one first peak very strong 
and stable at around 1 MeV, originated by a single transition, as we can see in the next figures. There is also a second 
peak at higher energy in Pb. This peak is more fragmented than the lower energy peak and moves to lower energy 
loosing strength as we move downward in the figure. In the prolate case both peaks are originated from transitions 
between negative parity states, reminiscent of the transition wh-n/2 — > vh§/2 in the spherical case. 

Now, we can particularize the discussion case by case analyzing in detail each isotope with regard to its deformation 
signatures. In the next figures (Figs. 6-11) we show for each isotope a comparison of the GT distributions obtained 
from the spherical (upper panel), oblate (middle panel) and prolate shapes (lower panel). We also show the discrete 
spectra, not folded, as they come directly from the QRPA calculations. 

In Fig. 6 for 184 Pb we can see a bump between 2.5 and 3 MeV in the spherical case. In the oblate case we get 
a three bump structure with peaks centered at 1.5, 3, and 5 MeV, respectively. The two peaks at lower energy arc 
fragmented within a range of 1 MeV while the third peak is little fragmented. On the contrary, in the prolate case 
we see a very strong single peak at 1 MeV and a second narrow peak at 3 MeV followed by a tail spread up to 5 
MeV. Similarly, we can see in Figs. 7-10 the following characteristics which are common to all the isotopes: i) an 
accumulation of the strength between 2.5 MeV and 3 MeV, depending on the isotope, for the spherical shape; ii) a 
very fragmented bump structure spread from 1.5 MeV up to 3.5 MeV, as well as a single peak located between 4 
MeV and 5 MeV, depending on the isotope, in the oblate case; iii) a single very strong peak at very low energy and 
a more spread structure beyond 3 MeV in the prolate case. Fig. 11 shows the results for 194 Pb, where only a stable 
solution, which is spherical, is found with Sk3 (see Fig. 1). The strength in this case is located beyond the Qec 
value and therefore of no significance from the (3 decay point of view. Our calculation with Sk3 predicts a stable 
isotope although the half- life has been measured to be T\/ 2 = 720 s. We shall come back to this point later on when 
discussing our results for half-lives. 

Now we focus our attention on the isotopes 188 490,i92p| 3 ^ j-, ecause several features combine to make them the best 
candidates for /3-decay strength distributions measurements [22]. These features are: i) a dominant f3 + /EC decay 
component (see Table 1); ii) a relatively large half-life (T"i/2 = 25 — 210 s); and iii) a large enough Qec energy (3.3-4.5 
MeV). In the case of 188 Pb, practically all the strength is found below its Qec value except for one peak at 4.5 MeV in 
the oblate case. A signature of deformation in this isotope would be the existence of GT strength below an excitation 
energy of 2 MeV, a region where the spherical shape does not generate any strength. If a very strong narrow peak 
is found in this low energy region, this would be an indication of a prolate shape, while if very fragmented strength 
is found this would indicate an oblate shape. The situation is similar in 190 Pb, but here it is even more evident the 
existence of three regions. Appearance of strength below 2.5 MeV would indicate again a deformed nuclear shape. If 
the strength were located below 1.5 MeV in the form of a narrow peak, it would indicate a prolate shape while its 
absence would rule out this possibility. In the case of 192 Pb, part of the strength is already located beyond the Qec 
window but still we get very clean signatures of deformation. The existence of a narrow strength distribution very 
close to the ground state would clearly indicate a prolate shape of the decaying nucleus and the absence of this strong 
peak together with the existence of strength extended up to 3 MeV indicates an oblate shape. The lack of strength 
below 2.5 MeV would rule out a deformed shape. 

We have also calculated the half-lives (Eq. (2.8)) corresponding to the (3 + /EC decays of these isotopes with the pur- 
pose of contrasting our results against the measured half-lives and checking that there is no fundamental disagreement 
with experiment. The f3 + / EC decay mode competes in this mass region with a-decay and the relative contribution 
of the two decays can be seen in Table 1 as the percentage of the (3 + / EC decay with respect to the total decay [28]. 
We can also see in Table 1 the experimental half-lives [28], as well as the calculated half-lives from the spherical, 
oblate and prolate shapes with the forces Sk3 and SG2, using quenching factors given in Eq. (2.9) and experimental 
Qec values. In principle, different Q- values should be used in the calculations for different starting/ending energy 
minima. Nevertheless, the Q-values calculated from the binding energies of parent and daughter are quite similar 
for the various shapes and close to the corresponding experimental values. We have also checked that the half-lives 
calculated with the experimental Q ec values change in less than a factor of two when Q ec is varied within 1 MeV 
around this value. This is at least the case in the more relevant isotopes 190 > 192 pb, where the [3 + percentage is above 
99%. One should also mention that the spherical cases are more sensitive to the Qec values because the location of 
the GT resonance close to the Qec value may be critical. Nevertheless, the half- lives in the spherical cases are always 
much larger than the corresponding to the deformed cases. 

The decay of the most unstable isotopes 184 > 186 pb are dominated by a-decay and therefore the predictions for the 
/3 + /EC decay half-lives should be above the measured value. This is indeed observed in the table. In the case of 184 Pb, 
the (3 + percentage is small (20(15)%). Then, we expect the calculated [3 + /EC half-lives to be roughly about one 
order of magnitude larger than the total experimental half-life and this is observed in the table. For 186 Pb the decay 
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is shared by a and (3 emission in comparable proportions. Thus, we expect (3 + /EC half-lives of the order of twice the 
total half- life measured, and this is roughly the case in Table 1. We also observe a general characteristic common to 
all the isotopes, the spherical shape produces always the largest half-life and the prolate shape the lowest one with 
the oblate half-life in between. The most interesting cases are lss > 190 > 192 p\). These decays are clearly dominated by 
/3 + /EC decay with half-lives ranging from T X j 2 = 25.5 s up to T\/ 2 = 210 s. In 188 Pb the experimental 7\/ 2 = 25.5 
s, is compatible with the prolate and oblate shapes in both Sk3 and SG2. In 190 Pb the experimental T1/2 = 71 
s, is much closer to the oblate calculation (T\ji = 83 s) and about a factor of two larger than the prolate half- life 
(T\/2 = 27 s) in Sk3. The SG2 results are in this case below experiment. The measured half-life of 192 Pb (Tj/2 = 210 
s) is well reproduced by the oblate shape, but it is a factor of 8 larger than the half- life of the prolate shape in Sk3. 
The results from SG2 are a factor of 2 (3) below experiment for the oblate (prolate) shape. In the case of 194 Pb, 
as already indicated, the Sk3 force predicts only one minimum, which is spherical and stable against /3-decay. The 
oblate shape corresponding to the force SG2 predicts a half-life that agrees nicely with experiment, while the half-life 
obtained with the spherical shape is too large. 

As a general comment we can mention that the f} + / EC half-lives predicted by the spherical shapes are in general 
much larger than the experimental half-lives. As a consequence, even though the spherical shape corresponds experi- 
mentally to the ground state of the nucleus, according to our calculations, the [3 + / EC half-lives are determined to a 
large extent by the f} + /EC decay of the oblate and prolate shape isomers. Within the present calculations, we would 
be unable to reproduce the half-lives if we would restrict ourselves to spherical shapes. 

IV. CONCLUSIONS 

We have studied the energy distribution of the Gamow- Teller strength in the neutron-deficient 184 4 86 ,i88,i90,i92,i94pj :) 
isotopes. The study has been done within a deformed pnQRPA formalism with spin-isospin ph and pp separable 
residual interactions. The quasiparticle mean field includes pairing correlations in BCS approximation and it is 
generated from a deformed HF approach with two-body Skyrme effective interactions. 

The parameters of the Skyrme forces are previously determined by fits to nuclear matter and ground state properties 
of spherical nuclei and the pairing gap parameters are extracted from the experimental masses. The equilibrium 
deformation is derived sclf-consistently within the HF procedure. The coupling strength of the ph residual interaction 
is fixed to reproduce the energy of the GT resonance in 208 Pb. The value obtained is close to the coupling strength 
derived consistently from the Skyrme force. Finally, the coupling strength of the pp residual interaction is taken from 
existing global parametrizations, although we find our results to be quite insensitive to this value. Our calculation 
has basically no free parameters to fit locally any of the particular isotopes discussed. 

Although the present procedure to generate the structure properties may be enough for a first estimate of the 
GT strength distributions, it would be interesting in the future to improve this description and discuss the effect of 
using the most recent Skyrme forces available in the literature together with a treatment of the pairing correlations 
using zero-range density dependent pairing interactions, where issues such as renormalization or regularization of the 
pairing force can be analyzed in this context. 

We have investigated the sensitivity of our results to the forces used and we arrive to the following conclusions: 
i) In a given isotope, the energy deformation curve and, in particular, the relative energies of the spherical, oblate, 
and prolate minima, are very sensitive to the Skyrme and pairing forces used. The ground state shape predicted 
may be modified with different choices for these forces. Nevertheless, the qualitative behavior of the energy profile in 
each isotope remains unchanged against changes in the Skyrme and pairing forces in the sense that the deformations 
at which the minima occur are hardly shifted and that the shape coexistence between spherical, oblate and prolate 
shapes is a constant feature in the more neutron deficient isotopes, ii) The GT strength distributions calculated at 
the equilibrium deformations show specific characteristics for each deformation that remain against changes of the 
Skyrme and pairing forces. The effect of the deformation on the GT strength distributions is much stronger than 
the effects coming from the Skyrme or pairing forces used, iii) As a consequence, we have been able to identify clear 
signatures of deformation on the GT strength distributions of the Pb isotopes studied. These signatures are related 
to the profiles of the GT strength distributions, which are peaked at different energies depending on the shape of the 
decaying nucleus. A detailed analysis of the best candidates to be measured, 188 . 190 > 192 pi ; ) ) has been performed and 
the characteristic features of their GT strength distributions have been identified. 

The calculated (3 + /EC half-lives show a strong dependence on the shape of the nucleus. The spherical shape 
produces always the largest half-life, which is far away from experiment. The prolate shape gives the lowest one while 
the oblate half- life lies in between and agrees better with experiment. The conclusion from our calculated half- lives 
is that in order to reproduce the experimental (3 + /EC half- lives, the decay of the spherical isomer is not enough and 
one has to consider the decay of the oblate and prolate shape isomers as well. 
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TABLE I. Percentage of the f3 + /EC involved in the total decay of the Pb isotopes, experimental [28] and theoretical half- lives 
[s], obtained with the forces Sk3 and SG2. Experimental Qec values have been used. 



A % (P+/EC) T*;* al cxp Tf (Sk3) Tf /2 (SG2) 









oblate 


prolate 


spherical 


oblate 


prolate 


spherical 


194 


100 


720 (30) 






stable 


642 




2053 


192 


99.9941 (7) 


210 (6) 


207 


28 


4385 


132 


84 


763 


190 


99.60 (4) 


71 (1) 


83 


27 


344 


53 


12 


195 


188 


90.7 (8) 


25.5 (1) 


32 


18 


85 


24 


22 


62 


186 


60 (8) 


4.82 (3) 


15.4 


7.7 


47 


8.2 


6.2 


16 


184 


20 (15) 


0.490 (25) 


6.4 


3.3 


22 


4.5 


3.2 


6.7 



9 




FIG. 1. HF energy of the isotopes I84,i86,i88,i90,i92,i94 pb obtained from constrained HF+BCS calculations with the force 
Sk3 and fixed pairing gap parameters as a function of the quadrupole deformation (3. 
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FIG. 2. HF energy of 184 Pb as a function of the quadmpole deformation f3, for various forces and pairing treatments. 
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FIG. 3. Half-lives for the oblate and prolate shapes of 190 Pb as a function of the coupling strength of the particle-particle 
residual interaction k pp . In the inset we can see the B(GT) strength as a function of the excitation energy for various shapes 
and Kpp parameters. 
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FIG. 4. Gamow- Teller strength distributions in Pb as a function of the excitation energy in the daughter nucleus. From 
left to right panels we show the results for oblate, spherical, and prolate shapes. In the upper panels we show the results from 
Sk3 (upward) and SG2 (downward) with fixed gap parameters, while in the lower panels we show the results from fixed pairing 
gaps (upward) and constant pairing strengths (downward) for Sk3. 
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FIG. 5. Gamow-Teller strength distributions in the I84,i86,i88,i90,i92 pb isotopes for sp h e rical (left), oblate (middle) and 
prolate (right) shapes. Results are obtained from Sk3 force with fixed gap parameters. 
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FIG. 6. Gamow- Teller strength distributions (discrete and folded) in 184 Pb. The experimental Qec energy is shown with a 
dashed vertical line. 
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FIG. 7. Same as in Fig. 6 for 186 Pb. 
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FIG. 9. Same as in Fig. 6 for 190 Pb. 
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FIG. 10. Same as in Fig. 6 for 192 Pb. 
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FIG. 11. Same as in Fig. 6 for 194 Pb. 
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